Renormalisation of flows on the 
multidimensional torus close to a KT 
frequency vector 



Joao Lopes Dias* 
Department of Applied Mathematics and Theoretical Physics, 
University of Cambridge 
Silver Street, Cambridge CB3 9EW, England 

April, 2001 



Abstract 

We use a renormalisation operator 1Z acting on a space of vector 
fields on T d , d > 2, to prove the existence of a submanifold of vector 
fields equivalent to constant. The result comes from the existence of 
a fixed point u> of 1Z which is hyperbolic. This is done for a certain 
class KTd of frequency vectors cu £ R d , called of Koch type. The 
transformation 1Z is constructed using a time rescaling, a linear change 
of basis plus a periodic non-linear map isotopic to the identity, which 
we derive by a "homotopy trick" . 
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1 Introduction 



In this paper we consider a renormalisation transformation 1Z that acts on 
a Banach space of analytic vector fields X on T d = M. d fZ d , d > 2. Its 
domain is an open ball around a non-zero constant vector field uj £ M. d 
(that generates a linear flow). The operator 1Z is a C 1 infinite-dimensional 
dynamical system with a hyperbolic fixed point at u. We show that the 
stable and unstable manifolds, given by the local dynamics of TZ around 
uj, yield different equivalence classes of vector fields. The elements of the 
codimension-((f — I) stable manifold are equivalent to u, as TZ is designed 
to be an equivalence between flows. In other words, the renormalisation 
asymptotically contracts locally any X into a (d — l)-parameter family of 
constant vector fields that includes uj. The map TZ basically consists of the 
composition of a time rescaling with a diffeomorphism of T d isotopic to a 
linear automorphism in GL(d,Z). 

This idea of renormalising vector fields is due to MacKay II, and we 
follow an approach inspired by the work of Koch || on renormalisation for 
(i-degrees of freedom analytic Hamiltonian systems. The latter applies to 
the problem of stability of invariant tori associated with a frequency vector 

UJ. 

The linear transformation is the main feature of the operator. It is a 
change of the basis of T d lifted to M. d , made to enlarge the region around the 
orbits of the linear flow. By iterating this, we look closer at chosen regions, 
but periodicity and the whole torus are kept at each stage. 

We also need to rescale the time as orbits take longer to cross the "new" 
torus. The topological characteristics of the trajectories are not affected by 
such transformation. 

The non-linear part of 7Z, isotopic to the identity, is obtained by a ho- 
motopy method to reduce the perturbation. We use a flow of coordinate 
changes (diffeomorphisms) to find one which fully eliminates some of the 
Fourier terms of X. Those are chosen to be the easiest to do so, since we 
keep the "resonant" terms associated with "small denominators" (as it ap- 
pears in the usual KAM theory). We are not concerned with trying to cancel 
all the perturbation terms because the linear transformation shifts resonant 
to "non-resonant" terms. Eventually all perturbation terms are eliminated 
while iterating TZ. The separation between resonant and non-resonant terms 
appears in ||. 

If a vector field is a fixed point of the above procedure, it means that 
its orbits exhibit self-similarity between the different scales. This is trivially 
deduced for the fixed point uj. 

The existence of the linear part, based on an arithmetical condition []|], 
depends on certain conditions imposed on uj. In particular, for d = 2, that 
corresponds to the set of vectors with a "quadratic irrational" slope. Specif- 
ically for the two-dimensional case, it is defined in a family of renormali- 
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sation iterative schemes allowing a full Lebesgue measure set of diophantine 
vectors. The methods involved therein are not easily generalisable to higher 
dimensions, it is required a suitable choice of a multidimensional continued 
fraction expansion algorithm. 

This paper is organised as follows. In Section |2| we recall the basic ingre- 
dients of equivalence of flows, and in Section || we introduce the renormali- 
sation idea through an example. We rigorously construct the operator 1Z in 
Sections [|, ||, ^ and ^, with proofs of the statements also given in Section ||. 
Finally, in Section ||| we present the main result using the spectral properties 
of the derivative at the fixed point, with further discussion in Section [U]. 

2 Equivalence of flows 

Consider a continuous vector field X on T d that generates the flow fa, 8 = 
X(6), 8 G T d , with lift $ t to the universal cover Mr. Choosing a norm II - II 
in W 1 , we define wx(9o) = lim^oo < ^ ) t(^o)/|| < l ) t(^o)|| to be the winding ratio 
of X for the orbit of 9q G T d , if the limit exists and lim^oo ||<&t(#o)|| = oo. 
Otherwise, if <3?t(#o) is bounded we put wx(8q) = 0. Also, if the limit does 
not exist or if ||$ t (0o)|| is unbounded but does not tend to infinity, the 
winding ratio is not defined. 

We say that two flows fa,i^t'- T d ^ T d are C r -equivalent if there is a C r - 
diffeomorphism h: T d — > T d taking orbits of fa onto those of ipt, preserving 
orientation. We are allowing to have h(fa(8)) = VV(/i(0),t) (^(8)), where t(8, •) 
is a homeomorphism of M for any 8 £ T d . This is the same to say that two 
vector fields X, Y on T d are equivalent if t'X ah = DhY where h is a 
flow equivalence as above and r' the time derivative of r. We emphasise 
that this relaxation of the usual requirement that t be preserved (ho fa = 
ipt o h) provides more satisfactory equivalence classes for flows, since we are 
in fact mainly interested in qualitative (topological) properties of the flow. 
Every C r -equivalence is isotopic to a map with a lift to M rf in the group of 
the linear automorphisms of the lattice Z d with determinant ±1, GL(d, Z). 
The isotopy is given by periodic homeomorphisms on T d . These kind of 
coordinate changes preserve the d-torus structure and volume. 

The set of winding ratios of a flow on T d generated by X is wx = 
{wx(8): 8 G T d }, which is called the winding set Any automorphism 
T: R d -> R d induces the map f : S^ 1 -> S^ 1 given by x Tx/\\Tx\\, 
where S d_1 = {x eR d : \\x\\ = 1}. Let T be in GL(d,Z). Then the winding 
set is invariant up to the action of T, with a C r -equivalence h isotopic to a 
map with lift T. That is, if X' = {Dhy 1 X o h, then w x >(9') = Tw x (8), 
where 8' = h~ 1 (8). In particular, the winding set is preserved by transfor- 
mations isotopic to the identity. 
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Figure 1: The action of T and U on an orbit of the lifted flow of X. 

3 Motivating example for d = 2 

We start by motivating the renormalisation procedure with a simple exam- 
ple. Consider the case d = 2 and the linear flow described by the differential 
equation 9 = u, with u £ M? and 6 £ T 2 . Given an initial condition 
#o £ T 2 for t = 0, the solution of the flow is 4>t(6o) = 6q + ujt mod 1, t > 0. 
The motion is simply a rotation with frequency vector uj and winding ratio 

= cj/\\cj\\. If, for all non-zero integer vectors k, k ■ uj := fciWi + ^2^2 / 
(i.e. the slope UJ2/UJ1 is irrational, assuming u\ does not vanish), then all 
the orbits are dense in the torus (the flow is minimal). Otherwise, they are 
closed curves (periodic orbits). 

Considering a perturbation of a constant "irrational" vector field uj lifted 
to M 2 , we want to determine under which conditions there is still equivalence 
to uj. Let X be a vector field close to uj arising from a time- independent 
perturbation. We choose e.g. uj = (1,7) where 7 = 1+ 2 ^ is the golden 
ratio. The main idea is to perform a change of basis, from the canonical 
base to {(0, 1), (1, 1)}, enlarging the region of M 2 around the orbits of uj 
(see Figure |l]). This is achieved by the linear transformation: 8' = T~ 1 6, 
where T = £ GL(2,Z). The eigenvalues and eigenvectors of T are 

given by: Tuj = juj and Tf2 = — ^Q, where £1 _L uj. We also rescale the 
time, t' = ^t, because the orbits take longer to cross the new torus. This 
transformation does not affect the unperturbed vector field uj, as it is given 
by -fDT- 1 o T(6') uj o T(9') = ^T~ l uj = uj. 

We should now consider a coordinate change h(9') = Ux T(9') with U 
satisfying = Id. The fundamental requirement is that it has to cancel 
the growth of the wiggles in the orbits of X, which are enlarged by the use 
of T (see orbit in Figure |l|). The idea is to perform a non-linear, close to 
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the identity periodic coordinate transformation Ux, in order to remove as 
many perturbation terms as possible. 

There are some terms which are more relevant - "resonant". To under- 
stand what they are and how they appear, consider the vector field X in 
the form of its Fourier decomposition: 

9 = ^ + 8^ fke 2mk - 6 

with f k G C 2 and e > "small". We can approximately solve the above 
system using the unperturbed solution plus an order e term. So, on the 
universal cover, 

9(t) = o + UJ t + ef o t + e y /* e 2 W tt.(flo-K^) +0 ( e 2) 

fcez 2 \{o} 

for t > 0. The resonant terms are the ones whose index k is almost perpen- 
dicular to OJ. 

We need to add some conditions to u in order to avoid having small 
denominators in the solution of the flow. This is done in the usual KAM- 
type proofs by imposing a diophantine condition on u, i.e. the denominators 
| A; • u\ in ([!]) admit a lower bound. We will show in the following that it is 
enough to find U such that it eliminates only "far from resonance terms". 
This is so because of the extra linear change of coordinates T described 
before, that is responsible for a "shift" of resonant into non-resonant modes. 
In some cases, by iterating this process the orbits can be straightened. 

We define U for any non-zero vector oj. Naturally, a restriction on to 
will appear again in the procedure with T (see also ||]), since this type of 
results fails for some vectors with irrational slope, in particular for some 
non-diophantine vectors O. 



4 Space of Analytic Vector Fields 

The following is valid for any ci-dimensional torus T d , d > 2, for a class of 
frequencies u to be defined. The vector fields considered are inside a ball 
around uj in some adequate space, and can be regarded as maps of M. d by 
lifting their domains. We make use of the analyticity to extend to the com- 
plex domain, so we deal with complex analytic vector fields. We construct 
the renormalisation operator 1Z and we look at the spectral properties of its 
derivative at the fixed point u, to relate to the local dynamics. 

The use of analytic function spaces as the domain of the renormalisation 
functional operator 1Z is justified by the usefulness of 1Z being C . Other- 
wise the "picture" of 1Z being a dynamical system with stable and unstable 
manifolds of the fixed points would vanish. The problem lies in the com- 
positions that appear in 1Z, thus reducing the degree of differentiability of 
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the image of the renormalisation operator. Considering C°° functions would 
also be a possibility, but complicate the technical parts of the method. 

Let r > and the domain lifted to a complex neighbourhood of K d : 
V(r) = {# G C d : ||Im0|| < , 

where || • || is the £i-norm on <C d . That is, \\z\\ = Ya=i with | • | the usual 
norm on C. We will also be using the inner product: z • z' = z i z l- 

An analytic function /: T>(r) — > C d , 2-7r-periodic in each variable is 
represented in Fourier series as 

m = E f^ 2mk - e , 

with coefficients fk G C d . The Banach spaces (Ad(r), || ■ ||r) and (A' d (r), \\ ■ \\' r ) 
are the subsets of the set of these functions such that the respective norms 

ll/||r = £ ||M|e r||fe|1 and ||/||; = £ (1 + 27r[|fc||) ||M|e r||fc|1 
fcez d kez d 

are finite. 

Let /3 > be a fixed value in the following. Consider the vector fields of 
the form X(-) = u + /(•), where u> G M. d and / G ^(z 9 )- Assume that X has 
no equilibrium points, X (6) 7^ for any 6 G T d . A condition like ||/||p < ||u;|| 
is enough to assure that there are no equilibria and, for any 6 G T d , there 
is no 9' G T d such that X(0')/\\X{9')\\ = -X{9)/\\X{9)\\. So, we could not 
have two symmetric winding ratios for such a vector field, since it would have 
implied having at least two symmetric normalised tangent vectors. The set 
of such vector fields on T 2 corresponds to the class of the Poincare flows 
(fixed-point-free with only one possible winding ratio, see e.g. Q). We 
study these systems in more detail in Q. 

5 Elimination of the Far from Resonance Terms 

We will be dealing with the vector fields considered above decomposed in 
Fourier series: X{9) = X^fcez d Xke 2mk ' e . As in 0, it is important to distin- 
guish different classes of terms according to their respective indices. 

Definition 5.1 For a fixed value of a > and ip G C d , we define the far 
from resonance terms with respect to ip to be the ones whose indices are in 

/;(^) = {fcez ,i : |^-jfe| ><7||fc||}. 
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Similarly, the resonant terms are = 7L d \l (J (ip). It is also useful to 

define the projections 1^ (VO an( i t(V0 f° r an y vector field X by 

fcei^WO 

and I = 1^ (VO + (VO ^ s ^ ne identity operator. 

The existence of a nonlinear change of coordinates £/, close to the iden- 
tity, that eliminates the far from resonance terms I~ (u) of a vector field X 
in a neighbourhood of a non-zero constant vector field, is given by 

Theorem 5.2 Let < p' < p, u £ C d \ {0} and < a < \\u\\. If X is in 
the open ball B C A' d (p) centred at u with radius 

p — p' 3 — y/6 a I 

4vr ' 6 ' 

£/iere is an invertible coordinate change U : T>(p') — > 2?(p) c/ose to i/ie iden- 
tity, satisfying 

[I- (lo)]{DU)- 1 X o U = 0, and C7 = Id i/ = 0. 

TTie map U : B ^ [1+ (uj)]Ad(p') given by X i— > (DC/) -1 X oU is analytic, 
and the derivative at every constant vector field u> + ip in B is equal to I+(u>). 
Moreover, 

\\U(X) -u\\p> < 2^1 + max 1^(3-^/6), 6(^+2)^1^) \\X - u\\' p . 

In Section |9| we include a proof of the above theorem using the "homotopy 
method". It is worth noting here that since we are requiring to eliminate 
only the far from resonance terms outside I^(oj), we avoid the problem in 
the proof of dealing with small denominators (see Section p^ ). 

6 Change of Basis 

All the coordinate changes of T d are isotopic to a linear transformation in 
GL(d, Z) which preserves the structure and volume of the d-torus. That can 
be seen as a map on the frequency vectors space, acting on a vector field by 
shifting terms in the space of indices Z rf and changing coefficients. 

Definition 6.1 For a fixed non-zero vector w£R' ! with rationally indepen- 
dent components, i.e. k-u ^ for any non-zero integer vector k, assume that 
lu = (l,u;2, •• • ,u!,i) where its components are real algebraic numbers (roots 
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of non-zero polynomials over Q). The vectors cu, c G R\{0}, are said to be 
of Koch type, i.e. to belong to the set KTd if the algebraic extension of Q 
by the numbers U2, ■ ■ ■ ,Wd (the smallest field containing U2, ■ ■ ■ ,u)d an d Q) 
is of degree d. 

An important property, the existence of a linear change of basis T cor- 
responding to the above class of vectors, follows from 

Lemma 6.2 (Koch ||] — Lemma 4.1) A vector u is in KTd if and only 
if there is an integral d x d matrix T , such that < < ■ ■ ■ < IA2I < 
1 < |Ai| and Tuj = \\uj, where Ai,...,Ad are simple eigenvalues ofT. 
Furthermore, some of the (integral) matrices satisfying these properties are 
in GL(d, Z). 



Remark 6.3 For d = 2, the numbers u>2 that produce a degree 2 algebraic 
extension of Q are the quadratic irrationals (roots of non-zero degree 2 
polynomials over Q). These are characterised by an eventually periodic 
continued fraction expansion Q. For d = 3, the base for the algebraic 
extension of Q by a and (3 is {1, a, a 2 }, where a has to be a cubic irrational 
(root of a non-zero degree 3 polynomial over Q). Then (3 = c\ + C2« + c^a 2 , 
d G Q. 

In what follows, a vector u S KT^ and a corresponding matrix T G 
GL(d,1i) are chosen according to Lemma ^2. Also, ui G M rf is chosen to be 



the (unstable) Ai-eigenvector of T such that u> - u = 1. Notice that KT^ is a 
subset of the diophantine vectors DC{j3) with = — 1 — In | Ai|/ In | A2 1 (Hi ■ 
Corollary 4.2), i.e. there is a constant C > that satisfies \oj-k\ > 
for every integer vector k 7^ 0. In addition, Ilj=i l-^il = 1) where Aj are the 
eigenvalues of T, with Ai G R since it is the only one outside the unit circle. 
Also, u ■ uj^ = 0, with ojv> being the eigenvector of T corresponding to A,, 
2 < j < d, because w • = T*oj ■ T~ l u^ = ■ We denote 

the transpose matrix of T by T*. 

From now on, we simplify the notation writing l£ and instead of the 
sets and the projections I^(uj), respectively. 

The next proposition determines that the change of basis corresponding 
to the matrix T above, acting on the space of "resonant vector fields", is 
analyticity improving. 

Proposition 6.4 Let < k < 1 and < p' < p. If up < p' , then, for some 
< a < ^||ci;|| _1 , any X G I^Ad(p') has an analytic extension on TT>{p). 
The linear map T(X) = T~ l XoT from I+Ad(f/) to A' d (p) is compact. 

The proof is given in Section [Oj. The values of a have to be chosen such 
that 

1+ C I K := {x G R d : \\T*x\\ < k\\x\\} 
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for a given k < 1. The proof of this proposition only requires that the 
resonant terms are inside the set I K . Thus, we could have different definitions 
for I + . The different choices are restricted by the part of the proof of 



Theorem 5^2 in Section that requires the denominator terms \lj ■ k\ to 
have at least a ||fc||-linear lower bound on I~ ■ 

It is also important to remark that one assumes that the orthogonal 
hypersurface in M. d with respect to u is contained in I K . This condition 
is indeed verified by choosing appropriate norms on M. d or by considering 
a sufficiently large power of T instead. The indices of resonant terms are 
therefore shifted towards I~ , eventually becoming far from resonance terms. 



7 Renormalisation Operator 

We are now in condition to construct a renormalisation operator 1Z based 
on the coordinate transformations introduced before. 



Fix u> G KTd, T and \j, 1 < j < d, as in Lemma |6.2| . We denote by 
K(X) the average on T d of a vector field X, i.e. 

E(X) = / X(6)d6, 
where d6 is the normalised Lebesgue measure on T d . 



Definition 7.1 Consider the maps U and T given by Theorem ^2 and 
Proposition |6.4| , respectively. The renormalisation operator 1Z is defined to 
be 

K{X)= _ 1 X, X = ToU(X), 
uj ■ E(X) 

for vector fields X. 

The rescaling of time is chosen to preserve the term u>, because without 
it the iterates of u, (T oU) n (u>) = A^^cj, n > 0, would go to zero. We could 
fix that by simply considering the product by Ai inside the calculation of 
X, but then all the elements of the one-parameter family X^ = (I + fi)uj, 
fj.EC, would be fixed points of the renormalisation. This "neutral" direction 
is contracted by the operator 1Z given in the definition, allowing uj to be an 
isolated fixed point. 

This choice of rescaling has the disadvantage of reducing the domain of 
the operator (see Proposition [7.2| below), when comparing with the size of 
the ball B obtained in Theorem |5.2| . In any case, this can be solved by using 
the rescaling suggested above, thus working with a one-parameter family of 
fixed points and a larger domain. 



Note also that, by Proposition 6.4, a < ^Mll ^ ^IMIj where 



1 2 represents the Euclidean norm. This means that the radius of an open 
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Figure 2: Steps of the renormalisation operator 1Z for d = 2, on the lattice 
Z 2 , the indices of the Fourier coefficients (notice that each coefficient is a 
vector in C d ). 



ball given by Theorem 5.2 can be taken to be Cct 2 ||u;|| _1 , assuming at least 
C < 4. Hence, the domain of 1Z only contains equilibria-free vector fields. 

The change of basis given by T turns some of the resonant terms into 
I~ where they will be completely eliminated by the operator IA of the next 
iteration of TZ. The use of such U is not essential, it would be enough to 
"sufficiently" reduce the resonant terms in a way that analyticity would not 
be lost. As the elimination is complete, we get a faster convergence to the 
fixed point and a much simpler and clearer analysis starting from the fact 
that the derivative of U is straightforward. If we would only apply U, all 
the neutral directions that are eliminated by the use of T would persist (no 
elimination of the resonant terms), and the analyticity domain would not 
be fully recovered. The right combination of the two coordinate changes is 
responsible for the usefulness of 71 (see Figure § for the d = 2 case). 

Proposition 7.2 Let p > 0. The renormalisation operator TZ is a well- 
defined analytic map from an open neighbourhood B C A' d (p) ofuo, to A' d (p). 

Proof: Consider the complex- valued continuous functional F{X) = Q ■ 
E(X), and restrict its domain to a neighbourhood B C A' d (p) of Aj -1 ^ such 
that F(B) C {z G C: \X±z — 1| < 1/2}, i.e. F is bounded away from zero in 
B. 

The domain B of 1Z is a subset of B (the domain of U) such that T o 
U{B) C B. That is satisfied for a small enough choice of the radius of 
B C B, since U(X) = u + l+f + C(||/||^ 2 ) and \XiF(T o U{X)) - 1| = 
|Ai(w • T^Ef) + C(||/ir/)|, with X = u + /. It is then sufficient to have 



|cD-E/ + 0(||/||^)| < \. Notice that Q -T^E / = T- U u>-Ef = A^(cj-E/). 



Hence, as F is bounded and analytic in B, the above, Theorem 5.2 and 



Proposition (Ll prove the claim. □ 
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8 Hyperbolicity of the Fixed Point uj 



The transformation 7Z was constructed in a way to hold an isolated fixed 
point at uj. The local behaviour of 1Z around it is given by the derivative, 
with which we hope to characterise the vector fields close enough to uj, 
realising the existence or not of an analytic equivalence. 

We can rewrite 7Z by the expression 1Z = T o T o IA, where the time- 
rescaling step is given by T(X) = [u>-E {X^X and DF^uj)/ = A^ x [/ — 
uj ■ E (/)] uj. The derivative of 1Z at uj is then the linear map 

DK{uj)f = £f-{uj-E(£f)]uj, feA' d ( P ), with Cf = X 1 ToI+f. 



This operator is compact by Proposition |6,4| and using the fact that 1+ is 
bounded. The eigenvalues of C are zero and those of T~ l multiplied by Ai. 
The vector uj gives the subspace corresponding to the eigenvalue 1. Since 
the renormalisation should not depend on the size of uj, this direction is 
fully eliminated in D1Z(uj) by the rescaling of time: D1Z(uj) uj = uj — (uj ■ 
uj) uj = 0. There are (d — 1) unstable directions of D1Z(uj) given by the other 
eigenvectors of T -1 , with eigenvalues XiXj 1 , 2 < j < d, of modulus greater 



than one (see Section [LC]) . 

The remaining eigenvalue of D1Z(uj) is zero, meaning that the elimination 
of only the far from resonant terms in conjuction with the change of basis 
are sufficient to eliminate all the non-constant Fourier terms of a vector 
field close to uj. Therefore, there is a codimension-(d — 1) manifold inside a 
neighbourhood of uj in A' d (p) being mapped by 1Z into itself, and we have 
proved 

Theorem 8.1 If ' uj € KT d , then the constant vector field uj is a hyperbolic 
fixed point ofTZ with a local codimension-(d — 1) stable manifold W s (uj) and 
a local (d — 1) -dimensional unstable manifold W u (u>). 



We have determined an analytic equivalence between the elements of 
W s (uj) and in particular with uj, that is given by the operator lim lZ n . 

n— ++oo 

The local stable manifold W s (uj) is the set of all X in some neighbourhood 



B with winding set wx = uj/\\uj\\ (see Section 10). On the other hand, 
the local unstable manifold W u (uj) is the affine space of the vector fields in 
the form X = uj + v £ B, where v G span-fa/ 2 - 1 , . . . ,uj^}, as this spectral 
subspace is 7£- invariant. A schematic representation of the renormalisation 
scheme described is in Figure S. 



9 Proofs of Theorem |5.2| and Proposition |t>.4 

9.1 Preliminaries 

We include here some technical details useful for the proofs. 
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W u {u>) 




Figure 3: The schematic renormalisation picture given by the action of 1Z 
in a neighbourhood of u. Note that dim W u (u) = d — 1 whilst W s (uj) is an 
infinite dimensional manifold (of codimension d — 1). 



Firstly, note that by writing / = (f 1 , . . . , f d ), one can rewrite the for- 
mulas of the norms on the spaces Ad{r) and A' d (r). In fact, 

d d 
11/11, = £11/% f GilM and ||/||; = £||r||;, ft A^r), 
i=i i=i 

where 

d 

ii/ii; = ii/iir-+En^ii- i = i,---,d. 

One has ||/|| r < ||/||;, A' d (r) C Mr) and \\f\\ r , < \\f\\ r , for < r' < r. 

Define the inclusion maps X: A d {r) — > -4,d(V) and X': .44(r) — ► «4' d (r'), 
< r' < r, by the restriction of the domain D(r) to D(r'), i.e. / i— ► /|x>(r')- 
It is easy to check that l{A d {r)) C A d (r') and J'(^(r)) C ^(V). The 
inclusion maps X and X' are compact linear operators. 

Consider a bounded linear map A: Ad(r) — > «4d(r) with (Af)(6) = 
A{6) f(6), where -A(0) is a linear map on C d . Its operator norm satisfies 
Pll < Ui =ill a ij'||r) where G -4i(r) such that, for each G £>(r), 
ciy(0) are the coefficients of the matrix which represents A{9) in the canon- 
ical basis. In particular, if A is the derivative of a function / G «4<i( r )> 
(£>/ 5 )(0) = Df(6)g(0), g G ^(r) and G 2>(r), then ||D/|| < ||/||' r . 

If a G ^li(r) and f,g£ A d (r), then a f £ A d (r), /■# G Ai(r), \\a f\\ r < 
||a|| r ||/||r and \\f ■ g\\ r < \\f\\r\\g\\r- Similarly, if a G A[(r) and f,g G A' d (r), 
then a/ G A' d (r), f-ge A'^r), \\a f\\' r < \\a\\' r \\f\\' r and \\f-g\\' r < \\f\\' r \\g\\'r- 
In this way, the spaces A d {r) and A' d (r) are Banach algebras. 

Lemma 9.1 Let < r' < r and f G A' d (r). If u G Ai( r ') IMIr' < 
(r — r')/4TT, then 

1. ||/(Id + «)|| r ,< H/ll 
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2. p/(id+«)ii<n/ir (r+r/)/2 , 

3. ||/(Id + «)-/[| r /<[|/[| / (r+r/)/2 H| r / J 
I \\Df(U + u)-Df\\< gj\\f\\' r \\u\\ r ,. 

Proof: Knowing that ||e 27rlfc ' u || r / < e 27r " fc """"'' / one gets [j] and The mean 
value theorem gives ||[ To prove |] apply again the mean value theorem now 
to \\djf l (Id + u) — djP\\ r > and also 

\\djf% + r')/2 = £ 2vrN(l + ||fc||)|/i|c r "l fc " e -(*- r ')ll fc "/ 2 

<^F £ (i + ii^ii)i^ie r||fc|| = ^7ii/ii;, 

fcez d 

where we have used the inequality: sup i>0 te~^ 1 < 1/(3 for /? > 0. □ 



9.2 Homotopy Method 

The sometimes called "homotopy trick" has been used in different problems, 
such as by Moser to show that all smooth volume forms on a compact 
orientable manifold are equivalent up to a diffeomorphism ([||] - Section 
5.1e). Other examples of its application are proofs of the Darboux theorem 
or the Poincare lemma (Q - Sections 5.5.9 and A. 3. 11, respectively), and 
Roussarie's proof of Morse's lemma. This procedure is also related to the 
"deformation method" used in different setups in KAM theory. It consists 
of a flow of symplectomorphisms that reduce in each iterate the size of the 
perturbation of an integrable Hamiltonian ]|, ||] . 

In the following fix a > 0, thus dropping the index a for the sets and 
corresponding projections P 1 . Choose also p and p' such that < p' < p. 
Let < 8 < \ and £ be a positive constant which will be determined along 
the proof and contains the restrictions on the size of the perturbation of uj 
depending on a and the norm of u, as will be seen in Section |9.3| . 

For vector fields in the form X = co + f, consider the open neighbourhood 
£ of the term /: 

£ = {feA' d ( P ): \\f\\' p <i}. 

The coordinate transformation U is written as U = Id + u, with Id as the 
identity transformation, and u in the open ball B of radius 5 in I~A' d (p'), 
i.e. 

B={uerA! d ( P '): \\u\\' pl <s}. 

For a given function / S £, consider the operator F ': B — > I~Ad(p'), 

F(u) = r (/ + Du)- 1 [u + /(Id + u)} , 

where / is the identity operator. This is simply the transformed vector field 
I' (DU)- 1 XoU, with X = uj + /. 
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Lemma 9.2 The derivative of F at u £ B is the linear map from I A' d (p') 
to l'A d (p'), 

DF{u) h= I-il + Du)- 1 [Df o (Id + u) h 

-Dh (I + Du)- 1 (u + f o (Id + «))] , 

with h chosen such that u + h G B. 

Proof: We need to compute the linear term on h of F(u + h) — F(u). As 
u £ B, we have that the bounded linear operator Du, from Ad(p') into itself, 
satisfies \\Du\\ < \\u\\' p , < 1. Hence, + = ^ Ii>0 (-Dji)™. Using the 

following formulas: 

n-l 

{Du + Dh) n = Du 11 + D u j DhDu 71 - 1 -* + 0(\\h\\ '/), n > 1, 

+oo n— 1 

£ £ {-\) n Dui DhDu n ~ x -i = -(I + Du)~ l Dh(I + Du)- 1 , 

n=l j=0 

and the Taylor expansion of / around + u(0), G £>(p), that gives 



/ o (Id + « + h) = f o (Id + u) + £>/ o (Id + u) h + 0{\\h\\', 2 ), 



we get 

+ 0O 

F(u + h)- F(u) = r J2 (-!)" [(Du + D/i) n (w + /) o (Id + u + h) 

n=0 

-Du n (uo + /) O (Id + u)\ 
= r (I + Du)- 1 [(uj + f) o (Id + u + h) 
-(w + /)o(Id + u) 

-Dh (/ + Dit)- ! (w + /) o (Id + u + /i)] 

+0(iwi;, 2 ) 

= I" (/ + Du)" 1 [Df o (Id + u)h 

-Dh (I + Du)- 1 ^ + f) o (Id + u)] 

That completes the proof. □ 
We want to find a solution for the equation 

F(u) = 0. (2) 

For that, consider a continuous one-parameter family of maps: U\ = Id + n^, 
A G [0, 1], with "initial" condition Uq = Id, i.e. no = 0, such that 

F(u\) = (l-A)F(uo). 
Differentiating the above equation in respect to A, we get 

DF{ux) ^dY = ~ m ' (3) 
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Remark 9.3 The derivative of F at u can be rewritten in the form 

DF(u) h = I~(DU)~ 1 [DX oUh-Dh (DU)~ l XoU]. 

Evaluating on voU, where v is the vector field that generates U, i.e. dll/dX = 
v o U and v o U = DU v, one gets 

DF(u)voU =r(DU)- 1 [v,X]oU 
= I-[v,(DU)- 1 X oU], 

where [v,w] = Dwv — Dvw is the commutator for vector fields v,w. Thus, 
writing X = (DU)^ 1 X o U, it suffices to solve the commutator equation 

r[v,x\ = -rx 

with respect to v satisfying l + v = 0, or, equivalently, invert DF(u). Note 
also that we allow to have v = v\. 

Proposition 9.4 Ifu £ B, then DF(u)~ l is a bounded linear operator from 
I~A d (p') to l~A' d (p') and 

WDFiu)- 1 ]] < ~. 

From the above proposition (to be proved in Section 9.3) we know that 
DF(u) is invertible for u E £>, thus we may integrate (||) with respect to A, 
obtaining the integral equation: 

u x = - f DFiu^Fitydix. (4) 
Jo 

In order to check that u\ £ B for any A £ [0, 1], we estimate its norm: 
\\u x \\' p ,< su V \\DF{v)- l F{Qi)\\' p > 

< supWDFiv^lWl-fWp, 
veB 

< T = *- 

Therefore, the solution of (^) exists in B and is given by (Q) when A = 1. 

Now, the open ball B as claimed is simply given by B = u + £ C A'(p) 
using e given by @. For X = uj + / G B we have \\u\\' p , = 0{\\f\\ p i). Thus, 

\\U(X)-u\\ pl < WY^i-DuYu + iDU^foUWp, 
n>l 




< 21 +m J^«(V0 + 2 
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Notice that U{X) = (DU)' 1 ^ + foU). Moreover, \\U(X) - uj - I + f\\ p > = 
This means that the derivative at uj of the map U is I + . Now, 
assume = 0([|/|| p /), u + V eBnC d , and write U{X) = U{lo + ^ + 

(f-m So, 



o 



p')> 



i.e. I + is the derivative of U at uj + tfj. That completes the proof of Theorem 
£2. 



9.3 Proof of Proposition |974 

To prove Proposition [O] we start by inverting the derivative of F at u = 0. 

Lemma 9.5 If f £ £ such that \\f\\' p < cr/4, i/te associated bounded linear 
operator I>F(0) _1 , /rom I~^4rf(/9 / ) to I~*4^(/c/), satisfies: 



° -4H/II' 



Proof: From Lemma |9.2| one has 



DF(0)h =I-(fh-Dhu) 
f-(D-u) 



h 



-I-f(D- uj)- 1 (D ■ u)h, 

where f h = Df h — Dh f. Thus, the inverse of this operator, if it exists, is 
given by 

DF(O)- 1 = —(D ■ lo)- 1 \l - rf (D ■ lo)- 1 ] . 
The inverse of {D ■ uj) is the linear map from l~Ad{p') to \~A' d {p'): 



(D ■ oj)- 1 g(0) = 



9k 



feei- 



So, 



HP- U )- 1 g\\L= £ 



2irik • uj 



l+27r||fc| 



2w\k-uj\ \\9k\\e- 



p'\\k\\ 



k£l 

^ Z_y 27TO- fc lly^H e 



p'IIA'II 



< 



1.91 



where the use of the definition of I (Definition |5.1|) was crucial to avoid 
dealing with arbitrarily small denominators. Hence, \\(D ■ uj)~ 1 \\ < — . 
Similarly, for /: l~ A' d {p') - A d (p'), 



||/%'<||/L 



+ 
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Hence it is a bounded operator with ||/|| < 2||/||' ,. Therefore, 



irfrD.ur'WK -\\ft <i, 



a 



assuming < a/A. So, I - I~f(D • lo)- 1 
in I~Ad(p') such that: 



-i 



is a bounded linear map 



I -!-/(£>■ a;)- 1 



< 



(7 



completing the proof of the Lemma. 

It remains now to estimate the variation of DF with u £ B. 



□ 



Lemma 9.6 Let < 5 < min{^, Given u £ B such that \\u\\' p , = 5 < 

5, the linear operator DF(u) — DF(0), mapping I~A' d (p') into I~Ad(p') for 
f £ E with 



'=£<£, is bounded and 



\DF(u) - DF(0)\\ < 



1-5 

Proof: Lemma |9.2| gives us that 



4vr 4-2<5\ 2-5, 
+ 7r )£ + ~ M 



p-p' 1-5 



1 - 5 



[DF{u) - DF(0)] h= T~(I + Du)~ l [Df o (Id + u) h - {I + Du)Df h 

-Dh (I + Du)- l (to + /) o (Id + u) 
+(L + Du)Dh(uj + f)} 
= + Du)- l {[Df o (Id + u) — Df — Du Df] h 

+DuDh{u + /) 

-Dh (I + Du)~ l [f o (Id + u) - f - Du (w + /)]} 
= l-(I + Du)- 1 {A + B + C}, 

where A, B and C are each of the respective three terms in the previous 
sum. Thus, ||(I + -Du) -1 !! — (1 ~~ ll^ n ll) ■ Using Lemma 9.1, one obtains: 



\M P ' < 



and \\B\\ p i < (\\u)\\ + 



Air 
p-p' 



p\\ u \\ P ' + 



i II M/ 
p> \\ u \\p> 



\h\\ 



p' 



' ,. Finally, noting that 



\\Du{L + Du)- 1 \\ < ||D«||.(1- \\Du\l)- 1 , 
and from Lemma |9.1| again, one gets the bound for the third term: 



lew, < T 



(p+p')/2 



\ u \\p' + \\ u \\p' (IMI + 



p') 



□ 
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To conclude the proof of Proposition 9.4, notice that, for u € B and 
/€£: 



WDFiu^W < (WDFiO)- 1 ^ 1 - \\DF(u) - DF(0)||) _1 



4?r i 4-26 
P-P' 1-5 



r 



The last inequality is true if we choose e and 5 to satisfy 



e < 5 



(7 



25 



2 (1 - 5f 



\U)\ 



1 + 25 + 



4tt 4-2(5 
+ 



1 - 5 \ P - ft 1-5 



(5) 



and 



'<3-3 



1 + 



a 



2 u) 



2 1 
< -. 



An effective value of i for each a can then be determined by the minimum 
value of the upper bounds given in Lemma 9J3 and (||). That can be obtained 
from a specific choice of 5 given by 



5 = min 



P-P 



' 3-V6 a | 11 



< 



2 2 



1 + 



a 



2\\uj\ 



4vr ' 6 

as long as c/||w|| < 1. Therefore, it is sufficient to impose 

. fa V6-2 ^ 



2 1 
< -. 



e : = mm 



-o5 



4' 12 

where the following inequalities were applied to 

25 



(6) 



' <Jl+ ° 



1-5 



2 to 



- , „ (S-V6)a a 

and — M < — < — • 

(1-5) 2 " " 2 2 



9.4 Analyticity Improvement 



We want to prove Proposition |6.4j , i.e. that X o T is analytic in £>(p) and 
has bounded derivative. Let a > such that 



max 



|T*x| 



: ||x|| ^ 0, \u ■ x\ < <x\\x\\ > < k 



Using X o T{9) = £ 



2mT*k-e 



\XoT\\ p < J2 \\X k \\e 



p\\T* 
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and for the derivative we have 

\\D(XoT)\\ p <2tt ||T*A;||e-^l fc ll||X fc ||e(^ + ^')ll fe lle' , 'll fe ll 
fce/+ 
< f 4X\\p', 

by choosing < (3 < p' — Kp, and using the relation sup t>0 ie^* < l/£ for 
f > 0. These bounds imply that \\T{X)\\' p < (1 + 2itk/ (3)\]T- l \\ \\X\\ p ,. 

Let r > p such that p' > rn. What was done above for T>(p) applies as 
well to V{r). Therefore, one can decompose T = To J', where J : l + A(p') — > 
A'(r) as before, and I: A'(r) — > »4'(p) is the inclusion map = X|x>( p ). 

Note that J" is bounded and X is compact, thus completing the proof. 



10 Spectral Properties of D1Z{uj) and Invariant Man- 
ifolds 

As a remains fixed, we continue to drop the index of the sets 1^ and the 
projections I ± . 

First we study the behaviour of the linear map C whose spectral prop- 
erties are closely related to the ones of DTZ(uj). Notice that EA'(p) = C d 
and (I — E )A'(p) are invariant subspaces of£as£oE = Eo£. 

Lemma 10.1 If go E DC{0), for any (3 > 0, and o~\\u)\\ < \, one finds 
constants a > and b, c > 1 such that 



\C n+1 (I-E)\\ < b\le- acn \\£ n (I-E) 



n > 



Proof: For every n > 0, we define l£ as the set of indices k that verify 
(T*) n k {0}. Hence, 

a\\(T*) n k\\ > \uj ■ (T*) n k\ = \T n uj ■ k\ = |Ai| n |w • k\, for k G /+. 

This inequality implies that 

|wife| < a\\{uj -k)uj\\+a (X[ l T*) n [k - (u-k)uj] 

Thus, assuming o"||a>|| < \ and knowing that k — (u ■ k) to is the component 
of k € 1^ on the spectral directions of T* corresponding to the eigenvalues 
Xj, j = 2, . . . , d, 



\u-k\< 



a 



1 — <t|M 



A2 

Ai 



|fe — (lo ■ k) u\ 



As to is a diophantine vector of degree > 0, there is a constant C > such 
that 

|wfc| > C , ||fc||- /3 - 1 , 
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Figure 4: The set C / for the two-dimensional case, 
and we get a lower bound for the norm of k 6 /+, 

n 

(7) 



,24-/3 >C }"^ 



Ai 
A 2 



2(7 

Let us now define the operator 1+ : A!{r) — ► ^4'(/o), where we have chosen 
r > p. This is a projection for the indices in together with an analytic 
inclusion. Making use of the bound ([?]), the operator norm of I + follows 
from 



ii,t/n; = E( 1 + 27r ii A; ii)ii^ii eri|fc|le_( " p)pl1 ^ 



,-ac" || jru/ 
lr i 



with a = (r- / 9)[C(2 CT )- 1 (l-cj||tj||)] 1 /(2+/3) > and c = |AiA 2 > 1. 

We can write £ n+1 with respect to C n , including the operator I+, in the 
form 

_ E)/ = £ n (I - E)I+£/, 

where C: A'(p) — > A'(r) is C followed by an analytic extension. That is, 
C = \\T o 1+ where T: I + A'(p) C I + A(p) — > A'{r) given as in Proposition 



3.4. The norm ||£|| < \\T\\ is determined by 



\\Tf\\' r <b\\T-^ '" ' 



p 1 



for some r > p and a constant 6 > 1 from the proof of Proposition 6.4 in 



Section 9.4. So, the claim follows from 



ll^+^i-E^ii^ii^i-^iip+iiiiriiii/ir,, 

and [|r -1 || < |Ai|. □ 

So, the direction given by the non-constant terms I — E is a stable 
eigenspace with eigenvalue zero. The spectrum of £ o E is simply the one 
of T _1 . That is, the eigenvalues are AiA" 1 , j = 2, . . . ,d, and 1 correspond- 
ing to the eigenvector uj. Hence, there is a "neutral" direction given by u 
whereas the remaining ones corresponding to u/ 2 ) , . . . , oj^ are all unstable. 
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The spectral properties of DTZiuo) are easily related to those of C because 

DK{u) = (I — PJE)£, 

where P^ is the projection on the subspace spanned by uj, i.e. P w /o = 
(u> ■ fo) uj. Note also that C commutes with the projection P^E because u is 
a eigenvector of T, thus of C. So, 

(I-P W E)£ = £(I-PJE). 

Therefore, the projection P M = (I — P W )E determines the linear space 
¥ u A'(p) spanned by {u/ 2 ), . . . uj^}, which is the unstable eigenspace of the 
linearised map at uj. In a similar way, P s = I — P" is the stable part. 

We are now in a position to define and determine the stable and unstable 
invariant local manifolds of 1Z at the fixed point uj. The stable one, W s (uj), 
corresponds to all the vector fields X in some neighbourhood B such that 
7Z n {X) G B, for any n > 0. Similarly, considering the inverse iterates of 1Z, 
we define the unstable local manifold W u (uj). 

If X G W s (uj), the winding set is equal to that of uj, since it is preserved 
by 1Z. This is in fact true because ui is an eigenvector of the matrix T G 
GL(d, Z), the linear part of the transformation induced by 1Z on the domain 
of X. The winding set is preserved up to the action of T: x i— > Tx/||Tx||, 
by iterating TZ. If the winding set of a vector X G B is wx = <-<VIMI, then 
it is also in W s (uj), otherwise would approximate the unstable direction and 
wx would be different. 

The set of constant vector fields is invariant under TZ, i.e. 7Z(C d H B) C 
C d , as can be seen from the way the non-linear coordinate change Li was 
constructed in Theorem This allows us to conclude that the invariant 
unstable local manifold is W u (uj) = F U (B). 

Lemma 10.2 We can find a ball B' C A' d (p) around uj such that, for any 
X G W s (uj) Pi B' , there are constants K > and 9 G]0, 1[ yielding 

\\1Z m {X) - lo\\' p < K9 m \\X - u>\\' p , m > 0. 

Moreover, there exists a > 0, c > 1 and N > satisfying 

||(I - E)K n (X)\\' p < e~ acn \\X - uj\\' p + 9 2n \\X - uj\\' p 2 , n>N. 

Proof: Fixing 0<t'<t<#<l and choosing an appropriate norm 
|| • || on A'(p) equivalent to || • ||' it is possible to have \\DTZ(u) f\\ < t'\\f\\, 
f G P s (A'(p)), since DlZ(uo) is compact. 

We can write X = uj + / + g where / = P S (X - uj) and g = F U (X - uj). 
By the stable manifold theorem, there is a constant A > for which 

|| 5 || <A||/ + <7|| 2 , and ||/|| < ||/ + g\\ + A\\f + ff || 2 . (8) 
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(9) 



The image of X is then X' = K(X) = lo + /' + g' with 

/' = DK{u) f, 
Sf = DK{u)g + 0(\\f + g 

The analyticity of 1Z and the formula of the second-order Taylor remaining 
imply that there exists constants B, C > satisfying 

ll/'ll < t'\\fl nm 

W\\<C\\g\\+B\\f + gf. (10) 

By combining this with (9) we obtain 

II/' + 9'\\ < t'\\f + g\\ + (At' + AC + B)\\f + g\\ 2 . 

So, there is a radius r > such that 

-uj\\ < t\\X -uj\\ 

for \\X — u\\ < r. Iterating this inequality m times we obtain 

\\R m {X) - uj\\' p < K9 m \\X - to\\' p , 

if B' = {X £ -A'dip) '■ \\X — uj\\ < r}, that is the first of our claims. 

From Lemma |10.1| we can find a > and c > 1 and a sequence of positive 
integers k n < en, e > 0, such that 

\\(I-E)[DK(Lu)] kn \\ < e~ ac ", 

whenever n is sufficiently large. Given such n and X G W s (co) fl -B', the first 
part of the proof showed that ||7?. ,n (X) — u\\ < 5 = t m r, with m = n — k n . 
Consider a map F : (/, g) ^ (/', g') defined by (|), and let /„ = P s (^ m (X) - 
uj) and c/o = P"(^ m (X) - u). The pair (/ ,ffo) satisfies 

II /oil < 2<5, 
liffoll < AS 2 . 

Applying the inequalities (0), if 1 < A < 1/5, then ||pi|| < (C + 9B)AS 2 , 
using the notation (fi,gi) = F l (fo,go), i > 0. The fact that ||/j|| is bounded 
by at least 25 allows us to iterate F for k n times, obtaining: 

f kn = [DTZ(oj)] kn f , 
\\g k J < (C + 9B) k -A5 2 , 

as long as (C + 9B) kn A5 < 1. This is in fact verified if we choose e > such 
that (C + 9B) kn t 2m < 9 2n , or, more strongly, 

[t- 2 {C + 9B)] e < {6/tf. 

Now, K n {X) =U j + f kn + g kn and (I - E)K n (X) = (I - E)(f kn + g kn ). 
Hence, 

||(I - E)ft n (X)|| < 2e~ ac "t m r + A9 2n r 2 , 
which completes the proof, if n is chosen sufficiently large. □ 
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